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ABSTRACT 

We present an expression for the non-linear evolution of the cosmological power spectrum 
based on Lagrangian trajectories. This is simplified using the Zel’dovich approximation to 
trace particle displacements, assuming Gaussian initial conditions. The model is found to 
exhibit the transfer of power from large to small scales expected in self-gravitating fields. 
Some exact solutions are found for power-law initial spectra. 

We have extended this analysis into redshift space and found a solution for the non-linear, 
anisotropic redshift-space power spectrum in the limit of plane-parallel redshift distortions. 
The quadrupole-to-monopole ratio is calculated for the case of power-law initial spectra. We 
find that the shape of this ratio depends on the shape of the initial spectrum, but when scaled to 
linear theory depends only weakly on the redshift-space distortion parameter, 0. The point of 
zero-crossing of the quadrupole, k 0 , is found to obey a simple scaling relation and we calculate 
this scale in the ZePdovich approximation. 

This model is found to be in good agreement with a series of A^-body simulations on scales 
down to the zero-crossing of the quadrupole, although the wavenumber at zero-crossing is 
underestimated. These results are applied to the quadrupole-to-monopole ratio found in the 
merged QDOT plus 1. 2- Jy IRAS redshift survey. Using a likelihood technique we have 
estimated that the distortion parameter is constrained to be 0 > 0.5 at the 95 per cent level. Our 
results are fairly insensitive to the local primordial spectral slope, but the likelihood analysis 
suggests n - — 2 in the translinear regime. The zero-crossing scale of the quadrupole is 
k 0 = 0.5±0.1/iMpc'’ 1 and from this we infer that the amplitude of clustering is 

= 0.7 ± 0.05. 

We suggest that the success of this model is due to non-linear redshift-space effects arising 
from infall on to caustics and is not dominated by virialized cluster cores. The latter should 
start to dominate on scales below the zero-crossing of the quadrupole, where our model breaks 
down. 

Key words: cosmology: theory - large-scale structure of Universe. 


1 INTRODUCTION 

The Newtonian analysis of linear growth of perturbations in an 
expanding universe is a well -understood problem (Peebles 1980; 
Efstathiou 1990). The extension of this into the non-linear regime 
has proven more difficult, owing to the strong mode coupling that 
irises in gravitational collapse, and most progress has been made 
through the use of tf-body simulations. 

Actual observations of galaxies in redshift space are further 
complicated (and made more interesting) by redshift distortions, 
caused by peculiar velocities of galaxies along the line of sight. Again, 
the linear problem is relatively well understood. Kaiser (1987) has 
shown that linear redshift distortions take their simplest form when 
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expressed in Fourier space, at least if structure is far from the observer 
so that the distortions are essentially plane -parallel. Here the redshift- 
space Fourier modes are related to the real-space ones by 

6\k> = (i + foim) (i) 


where a superscript s denotes a redshift-space quantity, p k is the 
cosine of the angle between the wavevector k and the line of sight, 
and £, the redshift distortion parameter, is the dimensionless growth 
rate of growing modes in linear theory', which is related to the 
cosmological density fi by (Peebles 19S0) 



( 2 ) 


in the standard pressureless Friedmann cosmology with mass- 
to-light bias b. It is through measuring the distortion parameter 



to analyse the quadrupole distortion measured in the QDOT pj us 
1.2-Jy redshift survey. We summarize our conclusions in Section 5. 
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& that one hopes to measure the cosmological density parameter, 

n. 

Kaiser’s formula (1) is valid only in the linear limit, and for 
plane-parallel distortions, neither of which approximations is well 
satisfied in reality. Several authors have now addressed the problem 
of generalizing Kaiser's formula to the case of radial distortions, 
while retaining the assumption of linearity (Fisher, Scharf & Lahav 
1994; Heavens &. Taylor 1995; Fisher et al. 1995; Zaroubi & 
Hoffman 1996; Ballinger, Heavens &. Taylor 1995; Hamilton & 
Culhane 1996). 

The main purpose of the present paper is to extend the analysis of 
redshift distortions into the non-linear regime, retaining the plane- 
parallel approximation for simplicity. Our approach to the problem 
is motivated by the consideration that the density in redshift space 
may appear highly non-linear even when the density in real space is 
only mildly non-linear. For example, a region that in real space is 
just turning around, a mildly non-linear condition, appears in 
redshift space as a caustic, a surface of infinite density, which is 
thoroughly non-linear. This leads us first to work in Lagrangian 
space (Section 2.1), and secondly to adopt the Zel dovich (1970) 
approximation (Sections 2.2 and 3). The Zel'dovich approximation 
is in effect linear theory expressed in Lagrangian space, inasmuch 
as it approximates the trajectories of particles as straight lines with 
(comoving) displacements growing according to linear theory. 

Our approach follows that of Taylor (1993), who studied the non- 
linear evolution of the power spectrum. Comparable approaches 
have been used by Bond & Couchman (1987) to evolve the galaxy 
angular correlation function, by Mann. Heavens & Peacock (1993) 
to evolve the real-space correlation function of clusters, and by 
Schneider & Bartelmann (1995) to evolve the real-space 
(unredshifted) power spectrum. Our approach differs from that of 
Hivon et al. (1995), who applied a perturbation expansion in 
Lagrangian space to second order to calculate the redshift-space 
skewness. 

For simplicity, we assume throughout this paper that the density 
field is unbiased, 6=1. Generally, the effect of evolution, if 
continuity is assumed (which is a fundamental assumption of this 
paper), is to tend to drive the bias factor towards unity. Continuity 
implies that the ratio of galaxy to matter density (1 + 6)/(l + $ M ) 
remains constant in Lagrangian elements, and if (somehow) a linear 
bias 6 = 65 m is established at some early time when 6 and $ M are 
both small, then the ratio of galaxy to matter density must be close 
to unity, (1 + W M V(1 + 6 m )» 1. It follows that the bias will be 
close to unity at later, non-linear epochs when d M is no longer small. 
Conversely, if bias is in fact important at the present, non-linear 
epoch (as may well be the case), then it must be that the assumption 
of continuity must break down in the not too distant past. Indeed, it 
may be that continuity is violated on an ongoing basis. In the present 
paper we choose to ignore this thorny problem, and simply assume 
an unbiased density field. 

We begin in Section 2 by deriving equations which relate the 
Lagrangian and Eulerian descriptions of density, and we calculate a 
general expression for the evolution of the power spectrum. We 
invoke the Zd’dovich approximation and Gaussian initial condi- 
tions and consider some of the general features of the resulting non- 
linear power spectrum. Some exact solutions for initially power- 
law spectra are derived. In Sea ion 3 we calculate the power 
spectrum in redshift space, again in the Zel’dovich approximation, 
and obtain a number of analytic and numerical results for the 
observationally interesting ratio of quadrupole-to-mooopole 
power. In Section 4, we compare the predictions of the Zel’dovich 
approximation with S -body simulations, and we apply our findings 


2 REAL-SPACE POWER SPECTRA 


2.1 Clustering in the Lagrangian frame 


In the Lagrangian theory of fluid mechanics, the central variable is 
the displacement vector field, rather than the density field of 
Eulerian space. This displacement field, t(q.t), maps particles 
from their initial Lagrangian coordinate, q , to the Eulerian-space 
coordinate, x. via the relation 1 


*(qj) = q + i(q.t) 


(3) 


where r) is the integral of the velocity field along the world-line 
of the particle: 


i(q,i) = J Ai v[xiq.t'),t] . 


(4) 


The advantage of the Lagrangian formalism for gravitational 
collapse was first pointed out in the seminal work by ZeFdovich 
(1970). Because the coordinate system of Lagrangian space is non- 
linear, moving with the particles themselves, this approach takes 
the analysis of perturbation growth of the density field into the non- 
linear regime while conserving mass density. 

However, the most useful observable quantity in cosmology is 
the Eulerian density field, p(x), inferred either from variations in the 
galaxy distribution, or indirectly from distortions in the Hubble 
flow, from gravitational lensing, or from the microwave back- 
ground. Continuity, along with the assumption of uniform initial 
density, implies that the relation between the overdensity 
6(x) ■ [p(x) - p]/p in Eulerian space and the Lagrangian displace- 
ment field {($) is 

6(j) = -{(?)] - 1, (5) 

where 6 d(jt) is the Dirac delta function. The Fourier transform 
$(*) = Jd\e* x 6(x) (6) 

of the overdensity 6(x), equation (5), is 

**) = jdV* V* *-l). (7) 

which relates the Fourier space Eulerian density field 6<k) to the 
Lagrangian displacement field £{q). 

The correlation function of Fourier modes defines the power 
Spectrum P(k): 

M)6\l t')> = (2vfP(k)6 D (k - k ') , (S) 


where the angular brackets denote ensemble averaging and the 
Dirac delta function arises from translational invariance. Trans- 
forming to differential and centre-of-mass coordinates, we find that 
as a result of translational invariance only the differential terms 
survive. We can then express the power spectrum by the integral 
equation 

/>(*) = jdV a «((e a ^ (9) 

where the term {e Ur ' a ') is identified as the generating function of the 
differential displacement vector field, ), of points separated by 


1 Throughout this paper we take the coordinates x and q to be comoving 
coordinates, defined relative to the general Hubble expansion of the 
Universe. 
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distance q = q\ - q: - 

AW)* >-*<*>• OO) 

Equation (9) for the power spectrum is valid at all times, not merely 
in the linear regime. However, in order to solve this completely one 
needs to have an expression for the displacement generating 
function and its evolution. In practice we shall have to approximate 
this. In the next section we use equation (9) to calculate the 
evolution of the power spectrum in the Zerdovich approximation, 
under the assumption of random Gaussian initial conditions. 


2 J2 Evolution of the power spectrum in the Zerdovich 
approximation 

A careful analysis of the Zel’dovich approximation has been given 
by Hui & Bertschinger (1995), who consider it in the context of 
local theories of gravity. For our purposes we only need the result 
that the particle displacement field scales according to linear theory : 

f(qj) = D(r)t L (q) 


i D{t) [ 
~(27T) 3 J <: 


where f L (^ I s J he linear displacement field, and 6 L (k) the linear 
density field, defined at some suitably early time. This approxima- 
tion is then a first-order theory in the displacement vector field, 
extrapolated to arbitrary later times. By construction the Lagran- 
gian coordinate system preserves mass, but in the Zel’dovich 
approximation at the expense of not satisfying the Euler equation. 
In the Lagrangian frame the density field is usually evolved along 
fluid streamlines, flfafx), f], but equation (5) allows one to calculate 
the density in Eulerian coordinates. 

The assumption of a Gaussian initial fluctuation field implies that 
the initial displacement field, { L , was Gaussian. The Zerdovich 
approximation (11) further implies that the displacement field {($) 
remains Gaussian for all time. Equation (9) for the power spectrum 
involves the Lagrangian generating function of the differential 
displacement field, which for a Gaussian vector field is 

= exp <?)])■ (12) 

where $ij(q)m (IjiqOtjiqz)) is the displacement correlation func- 
tion. In the Zerdovich approximation, the displacement correlation 
function grow s as the square of the growth factor D{t) in the non- 
linear as well as linear regimes: 

iMf.O = [Ddrftutq ) . (13) 

The linear correlation function ifojjtq) of displacements is related to 
the linear power spectrum P L (k) by 

tfuj ■ Uli(? i)tv<f:» 

(!4 > 

which at zero separation, q = 0, is 

*uj(0) = (t% = | k~'-P L (k), (15) 

where ({ 2 ) is the one-dimensional dispersion of displacements. 

Inserting the Gaussian generating function (12) into equation (9) 
for the power spectrum yields (Taylor 1993) 

P(k) = | d 3 * _ i ) . (16) 
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This equation (16). along with equations (13) and (14), provides a 
non-linear mapping between the initial and evolved pow er spectra, 
where the degree of evolution is controlled by the linear amplitude, 
D{i ), of perturbations. 

The unit term on the right of the integrand of equation (16), w f hich 
arises from subtraction of the mean density, yields a delta function 
at k = 0 on integration. The unit term can be ignored for k * 0, and 
we tacitly drop it hereafter. Nevertheless, it can help in numerical 
integrations to retain the unit term when k is small (linear regime). 

So far we have not assumed that the density field is isotropic, 
although we have assumed that it is homogeneous. In Section 3 w f e 
will apply the power spectrum equation (16) in redshift space, 
where the density field is not isotropic. In the remainder of the 
present section we assume that the density field is unredshifted and 
isotropic. 

For an isotropic density field, the displacement correlation 
function t^iq) resolves into irreducible components parallel and 
perpendicular to the pair separation q: 


iMf) = Hq%q, + - q,q ,) . 

(17) 

where 


Uq) = jikP L (k)^kq) - 2 ^ 

(18) 

and 



(19) 

For a general potential flow, as here, the parallel variance is 
related to the perpendicular variance by (Monin & Yaglom 

1975; Gorski 19SS) 

Cq 

(20) 

which can also be derived from equation (14) with the definitions 
(17). It is convenient to define the differential displacement covar- 
iances and by 

vM<7W.(0)-i L(q). 

(21) 

t-iq) ■ i_(,q) - iMg). 
w hich are related by 

(22) 

6q 

(23) 

The (isotropic) power spectrum (16) can then be written in terms of 
$+ and 

P(k) = J d 

(24) 

w here m = £ * q Integrating equation (24) over the azimuthal angle 
of q about the wavevector k is trivial, and further integration over n 
yields 

P(k) = Tie riiuf&q >F[kq. k 2 $_ ), 

(25) 

where F(A,B) is the function 



(26) 

w ith erf(z) the error function (Abramowitz & Stegun 

196S) 

erf(z) ■ [ d/e“ r . 

tt"- Jo 

(27) 

The error function erf(z) for complex i is available 

in programs 
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Figure 1. The evolution of the dimensionless amplitude 
in the Zel'dovich approximation for random 
Gaussian initial conditions. The solid lines represent the non-linear power 
spectrum and the dotted lines represent the linear initial spectrum. 

such as mathematica (Wolfram Research). Equation (25) reduces 
the evaluation of the power spectrum to a single integral over 
separation q. Excepting some simple cases (Section 2.4), this 
integral must be done numerically. 

The integrand in equation (25) oscillates rapidly when kq is large 
and real, which can make the integral difficult to evaluate (Schnei- 
der & Bartelmann 1995). One way to resolve the difficulty is to 
continue the analytically into the complex plane, and to shift 
the path of integration over q into the upper complex plane, which 
converts the oscillations into exponential decay. Some experimen- 
tation may be required to determine a good integration path, along 
which the integrand is well-behaved. For example, in the case of an 
initially Gaussian power spectrum, Section 2.3 bclow y we found it 
satisfactory to integrate first along the real axis to q - min(2£>, 1/*), 
and then to complete the integration along a straight line slanted 
upwards in the complex plane, at angle tt/ 6 from the real line. 
Variants on this path, with or without the initial segment along the 
real axis, and with other slant angles, worked in other cases. 

23 Example: a Gaussian power spectrum 

As an example, in Fig. 1 we show the Zel’dovich evolution of a 
power spectrum that is initially Gaussian, 

P L (k ) = 2w J exp(-* 2 ) . (28) 

Tie quantity plotted is the dimensionless amplitude A(i), the 
square root of the dimensionless power spectrum defined by (Pea- 
cock 1992) 

^ . <*> 
Fig. 1 shows clearly the two main features of non-linear growth, the 
amplification of small-scale power from the formation of caustics, 
and the restriction on the growth of intermediate-scale power 
imposed by the positive density constraint in voids. 

Fig. 1 also illustrates the limitations of the Zel'dovich approx- 
imation at later stages. While in reality the small-scale density 
contrast increases following the collapse of structures, in the 
Zel’dovich approximation structures expand back outwards after 


collapse. Thus in Fig. 1 the amplitude A(*) reaches a maximum of 
order unity as a function of time, and thereafter declines. In reality, 
the amplitude A(k) would be expected to increase monotonically 
with time. 

2.4 Solutions for initial power-law spectra 

For initially power-law spectra, Pi(k) * k n t convergent results are 
attained in the Zel’dovich approximation for power spectrum 
indices in the range — 3<n<— 1. For ns- 3, the differential 
acceleration between particle pairs, hence the differential displace- 
ment field, receives divergent contributions from large scales, 
while, for «*-l, the acceleration diverges on small scales. In 
the latter case, non-linear processes are expected to intervene so as 
to truncate the divergence, but we do not consider such truncations 
here. 

For initially power-law spectra the differential displacement 
covariances are power laws: a q * 1 ■ It is convenient here to 

normalize the growth factor, D(t ), so that 

1M9) = 0 V"*’ (30) 

(a normalization based, for example, on the power spectrum is less 
convenient because it leads to divergence in the displacement 
covariance for n = —3 and —1)* The relation (23) implies 

\M<?) = (-n-l )vM<?) ■ 

With the normalization (30), the initial linear power spectrum P i(k) 
is [note that P L {k) is defined without the factor so 
P{k) = D 2 P L (k) in the linear regime] 

/> L (*) = 4tt( 2 - n) sin[(— n — l)ir/2]r( 1 — n)k* (32) 

which is positive over the interval -3 < n < -1, going through zero 
at n = -3 and -1. 

In the Zel’dovich approximation, initially power-law power 
spectra evolve in a self-similar fashion: 

A 2 (*,Z» = A 2 (D : *" +3 ) . (33) 

The self-similar evolution of scale-free initial conditions in the case 
of a fiat, 0=1, universe is familiar from the BBGKY hierarchy 
(Peebles 1980). In the Zel’dovich approximation, the property of 
self-similarity extends also to non-critical universes. 

We now give exact solutions for the cases n = — 1,-2, and -3 in 
the Zel’dovich approximation. 

2 . 4.1 n = — 1 

In the limiting case n — -1, the displacement field has an inco- 
herent Gaussian nature. This arises because the incipient diver- 
gence of acceleration on small scales causes particles at finite 
separations to acquire essentially random velocities relative to 
each other. However, if particles in an initially uniform field 
move about randomly, the resulting field continues to remain 
uniform. Therefore the power spectrum should remain infinitesimal 
even when the displacement covariance grows to a finite value. 
Mathematically, the linear power spectrum (32) goes to zero in the 

limit n-*-l- , . 

In reality, non-linear evolution will in this case invalidate 
linear extrapolation of displacements assumed by the Zel’dovic 
approximation. Nevertheless, it is interesting to consider n-~ 1 
a limiting case. Defining the small positive quantity y by 

i (34) 

ra -n-1 , 

which tends to zero as n — — 1, one finds that for = D 4 ’ 


C 1996 RAS, MNRAS 282, 767- 778 












equation (30), the dimensionless power spectrum, A 2 , evolves as 

& : (k,D) = 3yD 2 k 7 e~ D k . (35) 

The power spectrum (35) is that of an n = -1 Gaussian field, 
modulated by a random Gaussian distribution of displacements 
with a one-dimensional one-point dispersion of D. 

Although the power spectrum is itself infinitesimal, the ratio 
fi/Po of quadrupole-to-monopole redshift power, which we shall 
consider in more detail later, is well-defined and finite in the limit 


2.4.2 n = -2 


For an initial power-law spectrum of index n = -2, the integral 
(24) can be done analytically. Here = vL = D 2 q , and the 
dimensionless power spectrum evolves as 


A 2 (k,D) 


16 D 2 k 

ir(l + 4 D*k 2 ) 2 


3 *D 2 k ‘ 

+ 4(1 + 4 D*k 2 ) v2 a * 


(36) 


In the linear regime A 2 * 16D 2 */tt, while in the highly non-linear 
regime A 2 « (Z) 2 *)” 3 . 


2.4.3 /i = -3 


In the limiting case n — —3, differential accelerations between 
pairs are dominated by tidal contributions generated on large scales, 
inducing motions which appear locally like anisotropic Hubble 
flows. The displacement covariances are \p + m D 2 q 2 , = 2yfv 

The integral (24) can again be done analytically, yielding 

■ ,37) 

Here the power spectrum remains a power law, P(k) * k ~ 3 , at all 
times. As argued by Schneider & Bartelmann (1995), the persis- 
tence of the k~ 2 form results from the dominance of caustics in the 
density distribution. The amplitude initially increases, from an 
exponentially tiny value, as more regions reach the point of 
collapse, but then reaches a maximum, at D 2 «= 1/18, and subse- 
quently declines as £>~ 3 as regions pass through the point of 
collapse and continue streaming outward. 

We note that equation (37) also gives the asymptotic form of the 
evolution of the power spectrum in the Zel’dovich approximation 
on small scales whenever the spectral index is less than —3 on small 
scales. This is true for example in the case of the Gaussian power 
spectrum Pi(k) * 2ir 2 e~ 4 considered in Section 2.3. In this case, 
rms differential displacements between pairs are proportional to 
their separation on small scales. Define 


* 

r- o (f 60rr 


f 


d kk 2 P L {k) . 


(38) 


Then, for large k (larger than the wavenumber at which the spectral 
index passes through n = -3), 


C -1AI20 

which j$ the same as equation (37), with a in place of D. 


(39) 


3 REDSHIFT-SPACE DISTORTIONS 

In Section 2 we considered some of the properties of the density 
field evolving in the Zel’dovich approximation. Here we turn our 
attention to the properties of such a density field when viewed in 
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redshift space. We assume the Zel’dovich approximation through- 
out this section; in Section 5 we will test the validity of the 
approximation with N-body simulations. 


3.1 Zel’dovich approximation 

In redshift space, the observed radial position of a galaxy is given by 
its radial velocity, which measures not only the uniform Hubble 
flow, but also the peculiar motion of the galaxy. Thus the position s 
of a galaxy in redshift space appears shifted from its true position x 
by its radial peculiar velocity 5 : 

5 = x + [x • v(x)]x . (40) 

In the ZeKdovich approximation, the peculiar velocity v is related 
to the displacement { at all times by the linear relation 
v[xiq),t] s= f (qj) « ftiiq.t), where 0, equation (2), is the dimen- 
sionless linear growth rate. As discussed in the introduction, we 
assume an unbiased density field, b = 1. The position s of a galaxy 
in redshift space then appears displaced from its initial Lagrangian 
position q (this initial position q is the same in both real and redshift 
space): 

* = ? + {'(?); (41) 

this displacement is given by a redshift displacement t*(q) which in 
the Zei’dovich approximation is 

{*(?) = I (?) + * t(?)J* . (42) 

Equation (42) shows that the displacement of a galaxy in redshift 
space is stretched by a factor (1 + ft) along the radial direction. The 
overdensity 6 s in redshift space is then (compare equation 5) 

6 s (S) = | dV fi D (j - q - f ’(9)] - 1 ■ (43) 

The radial character of the redshift distortion ft(x • f )x complicates 
evaluation of the Fourier transform of the overdensity. However, 
matters simplify if the structure being observed is far away from the 
observ er, so that the distortions are effectively plane-parallel. In the 
plane-parallel approximation, the redshift displacement becomes 

t s (q) = i(q) + ftlz-i(q))L (44) 

where the line-of-sight direction z is taken to be fixed. The Fourier 
transform of the redshift density, equation (43), is then 

6*(k) s JdVc^'^c***** - 1) . (45) 

t 

It is convenient to introduce the vector A\ which is the wavevector k 
stretched by (1 + ft) along the z-axis: 

K = * + 0tf**)z, (46) 

whose magnitude is 

Kzzka+Wnl + tful)'* (47) 

with fi k the cosine of the angle between the wavevector and the line 
of sight 

Hk = Z-ic. (48) 

In terms of the vector K , equation (45) can be rewritten 

m = — 1). (49) 

Just as the displacement f 5 of a galaxy in redshift space is 
stretched by a factor (1 -f ft) along the line-of-sight axis z, equation 

2 From here on we work in distance units of the velocity, where the Hubble 
parameter H 0 = 1. 
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(44), so also is the displacement correlation function V'y ■ (£i tj) * n 
redshift space stretched by (1 + (3) factors along each 2 -component 
of the tensor: 

t'j = ($i* + & -tZkWui&ji + 

If the unredshifted displacement field { is Gaussian (as is true for an 
initially Gaussian field evolved in the Zel dovich approximation, 
Section 2.2), then the redshifted displacement field f 5 is also 
Gaussian. 

The derivation of expression (16) for the evolved power spec- 
trum in Section 2.2 did not assume isotropy, and remains valid in 
redshift space. Thus the redshift power spectrum P\k ) is given by 
equation (16) with the displacement correlation function taken to be 
the redshift-space quantity 


P>(k) - J d*q e* '* { , *' 1 - 1 ) . (51) 

As in equation (16), the unit term in the integrand can be dropped 
for k # 0, although again it helps in numerical integrations to retain 
the term when k is small. Symmetry about the line of sight z implies 
that the redshift power P* (A. pi) is a function of the magnitude k of 
the wavevector k and of the cosine ji* of the angle between the 
wave vector and the line of sight. 

The relation (50) between the redshifted and unredshifted dis- 
placement correlation functions and \l/ tJ allows the combination 
kjkjf# in the integrand of equation (51) to be written conveniently in 
terms of the vector A, equation (46), as k t k } 4s\ } = Then 

equation (51) for the redshift power spectrum becomes 


< 52 > 

where now n = ft and are the (unredshifted) quantities 
given by equations (21) and (22). Equation (52) is the redshift-space 
counterpart of equation (24). Integrating equation (52) over the 
azimuthal angle of q about the vector K yields a Bessel function Jq. 

P’ik) = J 2irqfydM - 7 1 ) ■ < 53 > 


where s and c are 




(54) 


■(l-l*)"- 


i + Pui 


, (55) 

d+2 eni+e-vl) 1 ' 2 

Except for the interesting case 5 = 0, pursued further in Section 
3.4, equation (53) cannoi be reduced further to any simple analytic 
form. Moreover, the oscillatory character of the integrand makes it 
difficult to evaluate numerically. Further details on the evaluation 
of the integral (53) are given in Appendix A. In Section 3.5 we adopt 


a different approach to evaluating the redshift power spectrum, 
which is to smooth it. 


32 Redshift power spectra for initial power tows 

We now give analytic solutions for the evolution of the redshift 
power spectrum in the Zerdovich approximation for initially 
power-law spectra P L * k m t for the cases n = — 1 and n = —3. 

In Section 2.4 we saw that the n = -1 spectrum gave rise to an 
incoherent random Gaussian velocity dispersion, resulting in a 
simple damping term for the non-linear power in the Zel dovich 
approximation. Extending this to the redshift domain, we can see 
that the effect will be an additional damping term in the line-of- 


sight direction - an approximation previously suggested by 
Peacock (1992) to model the effects of virialized clusters. With 
the normalization discussed in Section 2.4, we find that the redshift 
power spectrum is 

A s2 (k,D) = 3vD 2 k 2 {l + . (56) 

This is the same as the Peacock smoothing with a 2 = £>'(20 + 0‘), 
where o v is the one-dimensional one-point velocity dispersion. 
Notice that we also recover the linear Kaiser boost factor, 
(1 + fob 2 - Hence we conclude that this form has some justification 
beyond its phenomenological foundations for scales where the 
spectral index is — 1, and arises as a special case in the model 
presented here. 

For n = -2, there appears to be no analytic solution for the 
redshift power spectrum, unlike the unredshifted case, equation 
(36). 

In the case n = -3, we find the analytic solution 
t7 1 
d ( *’ D) = (l2^W+2^F+?^ 

xexp 

which is the redshift counterpart of the unredshifted solution (37). 


1 + 2 + 3 — 2fl~/Xji 

12Z> : (1 + 2&zj + $ 2 nb 2 


33 Multipole analysis of the redshift power spectrum 

The effect of redshift distortion on the statistical properties of the 
density field can be best interpreted by a multipole expansion of the 
redshifted power spectrum. The anisotropic spectrum can be 
expanded in Legendre polynomials: 

me 

/*(*) = ][V ( ( k)T ( (Kk). (58 > 

(m0 

where T ( (n t ) » the (th Legendre polynomial and n k is again the 
cosine of the angle between the wavevector k and the line of sight. 
The multipole moments of the power are given by 

m) = j <W PWPt (#**)• < 59 > 

Hamilton (1992) suggested using the ratios of the multipoles as a 
way of measuring 0, expanding the correlation function from 
Kaiser’s linear theory. Cole, Fisher & Weinberg (1994) developed 
this to study the anisotropy of the power spectrum, w here non-linear 
effects could be avoided by limiting the analysis to small wave- 
numbers. The distortion parameter, 0, could then be estimated from 
the ratio of monopole and quadrupole power: 

*$(*) _ llll? (60) 

r 0 (k) r+l/i+IF ‘ 

In the linear regime each observable mode contributes to an 
independent estimate of /S. 

Non-linear effects modify the quadrupole-to-monopole ratio 
from Ihe linear value (60). In the following section we derive an 
approximate expression for this ratio in the Zerdovich approxima- 
tion, and in subsequent sections we present numerical results for the 
ratio in the Zel‘dovich approximation, in AMx>dy simulations, and 
in observations. 


3.4 An approximation to the quadrupole-to-monopole ratio 

In this section we derive an approximation, equation (64), which 
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relates the ratio P^Pq of the quadrupole and monopole redshift 
power spectra to the evolution of the unredshifted power spectrum 
in the Zel'dovich approximation. 

In equation (53) for t he redsh ift power spectrum, notice that the 
Bessel function J 0 (kqs \J\ — p}) in the integrand equals one when- 
ever j is zero, which happens in the limit fi — 0 (equation 54) and 
also in the case m* = 0 or 3 where the wavevector is perpendicular 
or parallel to the line of sight. Moreover in the non-linear limit 
where * the K~ terms in the exponential cause the integrand of 
equation (53) to decay rapidly as q increases and again the Bessel 
function is sensibly equal to one (as is also the factor e^ M ) 
wherever the integrand is non-negligible. In all these cases, the 
redshifted power, equation (53), reduces to the same form as its 
unredshifted counterpart, equation (24), but with in the latter 
replaced by (AV*) 2 ^, where Kfk = 1 + (equation 47) given 
that s = 0 and c = 1 (equation 55). In other words, the redshifted 
power equals the unredshifted power evolved forward by 1 + fol 
in the growth factor D(t): 

P , (k.D)-P[k.(l+enl)D] . (61) 

As it happens, this same equation (61) is also valid in the linear 
regime, according to the usual linear relations P(k,D) = D : P L (k) 
and P s (k,D) = (1 + 0 ,i]) 7 P(k.D). Thus approximation (61) is 
valid (i) in the limit fi — * 0, or (ii) if p k = 0, i.e. if the wavevector 
k is perpendicular to the line of sight, or (iii) if ^ k = l, j.e. if the 
wavevector k is parallel to the line of sight, or (iv) in the highly non- 
linear regime, or (v) in the linear regime. Evidently the approxima- 
tion is of some generality. 

Let us consider more closely the limiting case 5 — 0, for w hich 
the approximation (61) becomes exact (within the context of the 
ZeLdovich approximation). Expanding equation (61) as a Taylor 
series for small 5 gives 

P‘(k.D) = P(k,D) + frlD OS - 0) . (62) 

Here the redshift power spectrum is the sum of a constant term and a 
term proportional to p k% which implies that the redshift power 
spectrum is a sum of monopole and quadrupole terms, for small 5 
The ratio of quadrupole and monopole powers for small 5 is, from 
equation (62), 

P 2 (k) 45 a In P(k,D) a 

n<k ) = 3 ai„(Z>-> (/3 ~ 0) • (63) 

Amongst other things, equation (63) predicts that the quadrupole 
redshifted power /^(X) goes through zero where the unredshifted 
power P[k f D(t)} reaches its maximum as a function of time in the 
Zerdovich approximation. 

We now propose a generalization of the result (63) to arbitrary 5, 
which is consistent with and to some extent motivated by approx- 
imation (61), and which accords with the numerical results of 
Section 3,6 below, which are further supported by the A^-body 
simulations in Section 4.1. The numerical results indicate that, at 
least in the case of initially power-law power spectra, the quadru- 
polc-to-monopoie ratio P^JPq appears to satisfy a simple scaling 
law with fit such that the shape of the ratio is insensitive to 5, while 
its amplitude is proportional to the usual linear ratio (strictly, the 
numerical results are for smoothed power spectra, but we assume 
that the results hold also without smoothing). We further find that 
the zero-crossing of the quadrupole scales with fi such that 
£>“ - 1/(1 + fi). These empirical results, combined with equation 
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(63), suggest the general approximation 

Plik) 504 -ijfi 2 B\nP[kAl + fi) Xf2 D) 

W*I+§ 5 + J5 2 a In {D-) (64 > 

The approximation (64) relates the quadrupole-to-monopole ratio 
PJPo of the redshifted power spectrum to the evolution 
a In P/d In D 2 of the unredshifted power spectrum in the Zerdovich 
approximation. Although we have numerical support for the 
approximation (64) only in the case of initially power-law spectra, 
we suggest that it is likely to be a good approximation for arbitrary 
power spectra. 


3.5 The generating function of smoothed pow er spectrum 
multipoles 


So far we have obtained a number of analytic results for the redshift 
pow er spectrum in the Zel ’dovjch approximation, but in general it is 
necessary to resort to numerics. Unfortunately, numerical integra- 
tion of the redshift power spectrum, equation (53), while feasible, is 
unpleasant (see Appendix A). One way to sidestep the numerical 
difficulties is to smooth the power spectrum. This is not such a bad 
idea because it is necessary to measure a smoothed power spectrum 
from observations. That is, the Fourier modes measured in a 
catalogue are already convolved with the window function of the 
catalogue, and one may choose to smooth the power spectrum 
further to reduce error bars. In addition, smoothing makes defining 
quantities such as the zero-crossing of the quadrupole, discussed in 
Section 3.6, more robust to random errors. 

The smoothed (unredshifted) power spectrum is defined as 

/*(*) = (65) 


We choose to adopt a power law times Gaussian smoothing kernel 


W s (k,k) = — k-/k-) 

2i7Jt JV+3 r[(2// + 3)/2] ’ 


( 66 ) 


w hich is analytically convenient, and goes to a delta function in the 
limit N — In subsequent sections we use N = 1. 

For the redshift power spectrum, the smoothed harmonics can be 
defined similarly by 


h(h = (26 + 1) Jd% P>{k)? t (MtWnUt.1 ) . (67) 

The ratio PoikyP^ik) of smoothed quadrupole and monopole 
harmonics satisfies the usual equation (60) in the linear limit. 

The whole hierarchy of smoothed spectra, for various N and (, 
can be derived from the generating function 


C(a,b) 


1 6 y k P , (k)e~* l '~ bk ' l ‘‘ . 


(68) 


In terms of this generating function, the monopole harmonic of 
the redshift power spectrum, smoothed with the window H)y , is 


H(b = 


T f C(fl,6 = 0) •> 

L ^J«^-H2ir* lv+J r[(2A r + 3)/2] J 


(69) 


while the quadrupole redshifted power smoothed with the same 
window is 


H(h = 


3 fi _ 3±1 LiT" ( CM) ] 
2 |3fl db\ M [ dflJ a ,i-:\ 2 iri lv+3 r(( 2 A' + 3 )/ 2 J J ’ 


(70) 
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Figure 2. Plots of the ratio /? of smoothed quadrupole and monopole powers 
against for initially power-law spectra with indices, from lower 

curs e to upper, of n = - 1 , —1.5, -2. -2. 5 and -2.9. in each case curs es for 
0 * 0 and 1 are shown, which lie close to each other (and practically on top 
of each other for n = - 1 and —1.5) when R is scaled to the linear amplitude 
R l and the zero-crossing scale k v . For each value of the spectra) index. 3 *= 0 
gives the slightly higher scaled ratio. 


Substituting equation (51) for the redshift power spectrum into the 
generating function (6S) gives 

C(a.b) = jd-Vtfte’* (71) 

where V,j is the matrix 

*,/<?) = a5 ‘j + bi n + ^j <0 > ~ ■ 02) 

The A'-part of the integral (71) is the Fourier transform of a 
Gaussian, which solves in the usual way to give 

= (73) 

where is the determinant of 

|*| = <•;* + a){(^ + + + fl)|^,(l + 0) : + a + b) 


aa0 + 0 2 )-b]) (74) 

ith fji = z-q and tfrziq) given by equations (21) and (22), and 
^ the <j<?-component of the inverse of ^, r 

-1 4-„ 

« — i^.| ' + a 


+(l-M : )[* + (2/S + /3 : ) + fc + M : ^_/S : ]) . (75) 

Evaluation of the smoothed monopole and quadrupole spectra from 
equations (69), (70) and (73) still involves a double integral over jx 
and q. However, the integrand here is well-behaved (at least for 
small A\ such as N = 1 ) and depends only on elementary functions 
[besides $Jq)]. This compares with the unsmoothed case (Appen- 
dix A), where the integrand is ill-behaved, and evaluating the 
harmonics involves a double integral over an infinite sum involving 
special functions. 


3.6 Numerical results for quadrupole and monopole power 

From the generating function of multipole moments, equation (73), 
we have calculated numerically the smoothed quadrupole and 


monopole moments of the power spectrum for initially power-law 
spectra P L * k" with various values of the index n . and for various 0. 
The smoothing kernel is U\ . equation (66), with N - 1 . Fie. 2 shows 
the ratio R of the smoothed quadrupole and monopole powers 


t P;(A) 

P 0 (k) 


(76) 


for a representative sample of cases, n = — 1. -1.5, -2, —2.5 and 
— 2.9, and for 0 = 0 and 1 in each case. Values intermediate 
between these give curves intermediate to those plotted. The plotted 
ratio has been divided by the linear ratio, 
R l « (|0 4- = 0 : )/(l + ^0 + |0 : ), equation (60), and the wave- 


number has been scaled to the zero-crossing scale, k 0 . 

It is immediately apparent from Fig. 2 that the shape of the 
quadrupole-to-monopole ratio R depends on the spectral index n . 
but is insensitive to 0 at fixed //. at least for n 5: -2.5. The 
insensitivity to 0 is an interesting result, and we have previously 
invoked this insensitivity in proposing the approximation (64) to the 
(unsmoothed) ratio. As to the spectral index. Fig. 2 shows that, for 
n £ —2, the ratio R actually rises up above its linear value before 
turning over and going through zero. Physically, the linear squash- 
ing effect is first enhanced by caustics in redshift space as structures 
approach turnaround, and is then subsequently negated by ‘fingers- 
of-God' as structures collapse. Fig. 2 indicates that for n ^ —2 the 
caustic effect wins over the finger-of-God effect at translinear 
scales. However, as the enhancement peak scales as 
£ pk - 0.2 1 ' n ~'"k u . we find that R decreases for n ^ -2. for fixed 
0 . in the regime k ~ k 0 . 

A useful fitting function to the curves in Fig. 2 is 




(k/k v ) 


r ,0AiM+3t 


j. 


(77) 


which is accuraie to within 10 per cent over the range 
0.1 < k/k 0 < 1. and -3 < « < — 1. We shall use the fitting function 
(77) in fitting the observed ratio R measured in the QDOT plus 1 .2- 
Jy survey to 0 and n, in Section 4.2. 

The zero-crossing scale of the quadrupole power also contains 
useful information about 0 and the spectral index. The zero-cross- 
ing scale discussed below refers strictly to that of the smoothed 
quadrupole power, but the zero-crossing scale of the unsmoothed 
quadrupole is not greatly different. 

Suppose that the zero-crossing of the smoothed quadrupole 
power PWk) is observed to occur at wavevector k 0 . Define a 
corresponding real scale by q (i m a/£ n , w here a is a fitting constant, 
discussed below . Then w e find that the dimensionless amplitude of 
the displacement correlation function 3l the zero-crossing scale q$ 
is fitted tolerably well by the following fitting formula: 


jMgn) 

<75 1 + 0 ’ 


■ 


(7S) 


In the Zd’dovich approximation, formula (78) w ith the values (we 
adopt here the value of a that best fits the A'-body simulations 
reported in Section 4.1, since the same value works well also for 
ZeKdov ich) 

a = 2.3 . C = 0.19 (79) 


is accurate to better than 20 per cent for initially pow er-law spectra 
over the observationally interesting range of spectral indexes 
-2s/f s-1 and linear growth rate parameter 0*0*1. As a 
matter of interest, in the case n = -1 the exact analytic result (in 
the Zel'dovich approximation) is t+iqotiql = l/[a : U + 0)}. Phy- 
sically, it is not too surprising that the zero-crossing of the quadru- 
pole should occur when the rms displacement ^]l 2 (q) between pairs 
has reached an appreciable fraction of their separation q. 


€' 1996 RAS. MNRAS 282. 767-778 



Non-linear power spectra 775 


The amplitude of the power spectrum itself at the zero-crossing 
k 0 is a more complicated function of spectral index. During linear 
evolution, the dimensionless smoothed monopole redshifted power 
$ik) is related to the dimensionless displacement correlation 
function ^Jq)/q' at q = a/k by (compare equation 32: note that 
the factor r[(5-f mJ/ 2]/T(5/2) below comes from smoothing) 

&S(k) = ^1 j ^(2-/?)sin[f-/r“l)Tr/2]ni-rt) 

TtO/2) (a/k) 2 

Non-linear evolution, which is starting to become important at zero- 
crossing, will modify the power from the linearly evolved value 
(80). Setting aside this effect of non-linearity, one concludes from 
equation (SO), combined with the approximation (7S), that the 
amplitude of the monopole power at zero-crossing is sensitive to 
the spectral index ti. but relatively insensitive to 8 [the 1/(1 + 8) 
factor in equation (78) tending to cancel the 1 + + \ factor in 

equation (80)]. 

It is worth emphasizing this interesting if somewhat disappoint- 
ing result: one might have hoped that the amplitude of the observed 
monopole power at the zero-crossing of the quadrupolc would 
provide a measure of $ in the mildly non-linear regime, offering a 
check on the value of 8 measured from the quadrupole-to-monopole 
ratio in the linear regime, and perhaps a test of the bias b. However, 
the hope is not realized: the power at zero-crossing is sensitive 
mainly to the spectral index, not to 8- 


4 COMPARISON WITH SIMULATIONS AND 
OBSERVATIONS 


4.1 Comparison with simulations 


Given the approximate nature of the Zel'dovich mapping assumed 
in Sections 2 and 3. it is important to check the accuracy of these 
results with A'-body simulations, where non-linear effects can be 
accounted for on smaller scales. The simulations that we ran used 
the AP^M code of Couchman (199J). generalized to low Q by 
Peacock & Dodds ( 1994 ). on a 64 ? grid with 64 ' particles. 

We ran simulations with power-law initial spectra with indices 
« = —1, —1.5 and —2. both for fi = 1. and for a low-0 spatially fiat 
model with 0 = 0.3 and a vacuum-energy contribution of Q \ =0.7, 
corresponding to 0 = 0.5. In each case three simulations were 
computed with differing random Gaussian initial conditions, IS 
simulations altogether. We assumed no bias in the density field, and 
evolved all the models for 300 time-steps so as to erase the effects 
of using the Zel'dovich approximation in the initial setting up of 
the simulation. 

Power multipoles were calculated by assigning densities to the 
grid by the cloud-in-cell algorithm and fast Fourier transforming the 
cube. The shot-noise component of the power was subtracted and 
the resulting Fourier space compressed by averaging over 
azimuthal angles. The multipoles were then calculated by a least- 
squares fit to each mode: 


x : <*) = ^2 "‘ (k) 




Y,Wk)T<<n,) 


(81) 


The weighting scheme used was the inverse variance. 
h \(k) = 1 fo'ik.n,). calculated from the intrinsic cosmic variance. 
a(k) = \^2^~ik). for each mode. The resulting power multipoles 
were then smoothed using the kernel function U' v with A' = 1, 
equation (66). and an inverse variance weighting scheme. 





Figure 3. Plots of the ratio H of smoothed monopole and quadrupole powers 
scaled by the zero-crossing. k 0 . From top to bottom on scales k < k 0 are 
plotted the numerical results for the Zel’dovich evolution for each power 
spectrum (solid lines), for 0 = 1 withn = -1. -1.5 and -2 (top three lines) 
and 8 = 0.5 with n = -1. -1.5 and -2 (bottom three lines). Overlaid are 
the results from a set of A'-body simulations with the same initial parameters 
(points). Errors on the simulations are from cosmic variance. 

Generally, we found that the Zel'dovich approximation ceased to 
provide a good approximation to the evolution of monopole or 
quadrupole power at moderately non-linear epochs. This is to be 
expected, since, as discussed in Section 2.3, the expansion of 
structures following collapse which occurs in the Zel'dovich 
approximation causes power to reach a maximum and then decline, 
whereas in reality powder continues to increase monotonically. 

However, as shown in Fig. 3, the Zel 'dovich approximation does 
appear to reproduce remarkably well the amplitude and shape of the 
ratio R of the quadrupole and monopole pow ers seen in the A'-body 
simulations, for scales down to the zero-crossing of the quadrupole, 
k s k {) . Small deviations occur at small A\ but these can be attributed 
to the lack of contribution from smaller wavenumbers in the 
smoothing process of the simulations. Small deviations are also 
seen in the n = -2 spectra, at least some of which may arise from 
cosmic variance and the influence of the finite box si 2 e on the longest 
wavelengths in the simulations. However, the change in shape of the 
quadrupole-to-monopole ratio with spectral index predicted by the 
Zel'dovich approximation is clearly seen in the simulations. All the 
models fail in the regime k > k 0 , where w e expect virialized clusters 
to provide the strongest redshift distortions (Jackson 1972). 

The scaling relation suggested by equation (7S) for the amplitude 
of the displacement covariance at the zero-crossing of the quadru- 
ple is also found to hold for the simulations. However, the 
amplitude C predicted by the Zel'dovich approximation was not 
so accurate. Applying a least-squares fit of the model to the 
ensemble of IS simulations, we found 

a = 2.3 ±0.4, C = 0.54 ±0.07 (82) 

with a formal \ 2 = 15.6 for 16 degrees of freedom. In general, 
while the form of the scaling relation (78) holds, the Zel'dovich 
approximation underestimates the wavenumber at zero-crossing of 
the quadrupole by a factor of approximately 2 in the cases 
considered. This can be understood as arising from the outflow 
of particles which follows caustic formation in the Zel'dovich 
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Figure 4. Plot of the ratio R of smocv.hed monopoie and quadrupole powers 
as a function of wavenumber k for the merged QDOT plus 3 2-.lv redshift 
survey. The solid line is a ftnoa model w iih 8 = O S and n — -2.5. and ith 
io = 0.5/iMpc' 1 . 

approximation. This outflow exaggerates the ‘finger-of-God effect, 
causing the zero-crossing of the quadrupole to appear at larger scales 
in the Zeldovich approximation than is actually the case. 

In fitting to observed data in the next section, we use the .V-bodv 
simulations to calibrate the zero-crossing of the quadrupole. 

4.2 Application to the QDOT plus 1.2*Jy redshift survey 

We have applied the results obtained above to estimate the distor- 
tion parameter. $. and the local primordial spectral index, /r, from 
the amplitude and shape of the ratio R of smoothed quadrupole and 
monopoie pow ers observ ed in the merged QDOT plus 1 .2-Jy IRAS 
redshift survey. We further infer the variance al of counts in 
8/r'Mpc spheres from the observed scale of the zero*crossing of 
the quadrupole. The observed quadrupole-to-monopoie ratio, taken 
from Hamilton (1995). is shown in Fig. 4. Superimposed on the 
observations is a model with 8 = 0.8 and ti = —2.5. The zero- 
crossing is estimated to be at In = 0.5 ± 0. 1 h Mpc 

We use the fining function (77) to fit the amplitude, which 
depends on 0, and the shape, which depends on ii. of the observed 
quadrupole-to-monopole ratio R. We construct a likelihood 
function 

L = exp{—a,M,jO,) (83) 

where a, = (data - model) is the difference between the ith data 
point and the model prediction. Af y is the data covariance matrix 
whose diagonal terms comprise the uncertainty on each data point 
and whose off-diagonal terms are estimated assuming that the 
underlying data have independent Gaussian-distributed errors, 
and the major source of correlations comes from smoothing. As 
our model breaks down on scales smaller than the zero-crossing 
point, and the errors on the observed quadrupole diverge at low 
wavenumber, we restrict our analysis to wavemodes in the range 
0.06 < klUt Mpc* 1 ) < 0.5. 

Fig. 5 shows the marginalized likelihood in the <«.£) plane. As 
n — ► —5, & is forced to infinity by our fitting function. The like- 
lihood function selects this range as the better fit to the data, ow ing 
to the slight inflection in the quadrupole-lo-monopoic ratio around 



n 

Figure 5. Likelihood function L (3- data), of 3 and the local spectral index, n. 
from the quadrupole-to-monopole ratio of the merged QDOT plus 1 2-Jy 
redshift survey. Contours are spaced b\ 6 In £ = -0.5 with the inner contour 
delineating the 6S per cent and outer contour the 95 per cent confidence 
region. While the spectral index is not strongly constrained. 8 < 0.5 is ruled out 
at the 95 per cent level. This rules out f? < 0.5 if //t4S galaxies are unbiased 

k = 0.2/iMpc” 1 . However, if we consider only the 95 per cent 
confidence limit given by the outer contour, low values of 8 < 0.5 
are strongly excluded. In terms of the galaxy bias parameter, this 
implies that the density parameter Q > 0 3b J '\ The value of the 
primordial spectral index on these scales is not strongly con- 
strained. but is consistent with //- -2. 

We estimate the clustering amplitude a s from the scale of the 
zero-crossing of the quadrupole, using equation (7S) with the 
parameters (82) calibrated from the ,V -body simulations. In Fig 6 
we plot the marginalized likelihood function £(a h ) = Jd3d/r 
£<0.«.ff S !data). We deduce that = 0.7 ±0.0 5 where we quote 



Figure 6. Marginalized likelihood function £(o s ) of the variance in 8/t 1 Mpc 
spheres inferred from the zero-crossing of the QDOT plus 1.2-Jj quadrupole- 
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Given the indirect method of our measurement, wc find this very 
encouraging. 

5 SUMMARY 

The main aim of this paper has been to study the redshift distortion 
of the power spectrum in the moderately non-linear regime. Our 
approach is motivated by the consideration that structures may 
appear in redshift space to be more non-linear than they really are - 
for example, regions that are just turning around in real space 
appear as caustics in redshift space. This leads us first to work in 
Lacrangian space, and secondly to use the Zeldovich approxima- 
tion, which is in effect linear theory expressed in Lagrangian space. 

We started by deriving an expression relating the power spectrum 
to the Lagrangian displacement field. We used this to determine the 
evolution of the power spectrum in the Zel'dovich approximation, 
first in real space, then in redshift space in the plane-parallel 
approximation. We presented some analytic solutions for initially 
power-law spectra. In particular, we showed that a spectrum with 
index // = — 1 gives rise to an incoherent Gaussian displacement 
field, producing non-linear redshift distortions with the same form 
as the Kaiser-Peacock model. 

We derived various analytic and numerical results for the 
observationally interesting ratio R of quadrupole and monopole 
redshift powers, whose value in the linear regime, 

R l = + + |£ : ). yields a measure of the distor- 

tion parameter 0. In the Zel’dovich approximation, the amplitude of 
the ratio R is set by the linear value R L . but the shape of R as a 
function of wavenumber k depends mainly on the spectral index n, 
and is insensitive to 0 at fixed n. The zero-crossing of the quadru- 
pole power occurs at the point where the Zel'dovich power 
spectrum is a maximum as a function of time. 

We have tested the Zel'dovich results against A- -body simula- 
tions with initially power-law spectra, in both 0 = 1 and low-Q 
models. The simulations show that the Zel dovich approximation 
ceases to provide a good approximation to the power spectrum at 
moderately non-linear epochs. Remarkably, however, the Zel do* 
vich approximation predicts rather well the amplitude and shape of 
the quadrupole-to-monopole ratio R on scales down to the zero- 
crossing of the quadrupole. when R is scaled to the scale of the zero- 
crossing. The Zel’dovich approximation underestimates the wave- 
number at zero-crossing by a factor of about 2, although it predicts 
correctly the way in which the zero-crossing scales with 0 and n. 

We have applied these findings to estimate the distortion para- 
meter /3, the local spectral index n. and the v ariance a * of counts in 
8/j -1 Mpc spheres, from the quadrupole-to-monopolc ratio/? mea- 
sured in the merged QDOT plus 1 .2-Jy redshift survey. We find that 
the distortion parameter is constrained to & > 0.5 at the 95 per cent 
level. The spectral index is not well constrained, but is consistent 
with at translinear scales. The clustering amplitude cr 8 , 

inferred from the scale of the zero-crossing of R , is 
a a = 0.7 ± 0.05, consistent with other estimates. 

The success of the Zel'dovich model in describing the quadru- 
pole-to-monopolc ratio R suggests that departures from the linear 
value R\_ at translinear scales are caused mainly by intall on to 
dusters, not by vitalized cluster cores. 
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APPENDIX A 

This section shows how to evaluate the ‘difficult' integral (53) for 
the redshift power spectrum in the Zel'dovich approximation. The 
integral is a two-dimensional integral over n = A' • q and separation 
q. We choose to develop the integral over analytically as an 
infinite sum. both because the numerical integral over y. is itself 
unpleasant, and because the expression as an infinite sum allows the 
rj-imcgrand to he continued analytically into the complex ^-plane, 
which stabilizes the integration over q . 

For small fc/s, the Bessel function J 0 {kqs\' 1 - m : ) the 
integrand can be expanded as a power series in 1 - /T. Integration 
of the leading term (which is 1) in the series leads to an expression 
similar to (25), and higher terms in I - can be generated by 
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repealed differentiation of this expression. The result is the infinite 
sum 


F s (k) = Re 



-q~dqe~ 




*± 

n«0 


i-fikqs/l) 2 ” 

n'r 



Fikqc.Bh 


(Al) 


where F is the function defined by equation (26), and s and c are 
given by equations (54) and (55). The sum in equation (Al) can be 
truncated to a few terms if kqs is small. For large k (non-linear 
limit), this truncation effectively works for all q % since then the 
exponential factor C “ A causes the integrand to become 

negligible when kqs is not small. 

The expansion (Al) fails when kqs is large and at the same time 
A'*(^ + \j/_ ) is small. Here a series expansion in leads to the 
following infinite sum: 


F\k) — Re | 4-q'dqt 




x 


x 

z 

nsQ 


n! 



i[.4 ; + ( kqs 


(A2) 


Expression (A2) works when is small and kqs is not small, 
which is complementary to the domain of validity of equation (Al). 
A program such as ma thematic a helps considerably with the 


algebra involved in taking the repeated derivatives in both equa- 
tions (Al) and (A2). 

The integration of equations (Al) and (A2) over q must be done 
numerically. As in the unredshifted case, the oscillatory character 
of the integrand for large kq can be converted into exponential 
decay by shifting the path of integration of q into the upper complex- 
plane. One should be careful here that, although the real parts of the 
integrands of equations (Al) and (A2) agree for real q , their 
imaginary pans differ, so that their analytic continuations in the 
complex <y-plane differ. The integrands could be made identical by- 
adding their complex conjugate expressions, but this would intro- 
duce an undesirable exponentially increasing component in the 
upper complex plane. Therefore, if it is necessary to use both 
expressions (Al) and (A2) in integrating from q — 0 to x. then the 
transition from one expression to the other along the path of 
integration must occur on the real line. 

In the cases that we tried, an integration path similar to that 
described at the end of Section 2.2 proved satisfactory. Namely, 
integrate first some way along the real axis, then complete the 
integration along a straight line upwardly slanted in the complex 
plane. In the case K’^Aq = 1 Iks) 2 1. equation (Al) works every- 
where on the path of integration. In the case K l $„iq = 1 Iks) < 1, 
equation (Al) works for the first part of the integration, along the 
real line, but equation (A2) is required for the second part, a suitable 
path being vertically upwards in the complex plane. 

Computing the quadrupole and monopole harmonics of the 
redshift power spectrum requires a final numerical integration of 
F\k) over n k with for ( = 0 and 2, equation (59). 


© 1996 RAS. MNRAS 282, 767-778 


